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Abstract 

We investigate the time evolution of waves in evanescent media generated 
by a source within this medium and observed at some distance away from 
the location of the source. The aim is to find a velocity which describes a 
causal process and is thus, for a medium with relativistic dispersion, limited 
by the velocity of light. The wave function consists of a broad frequency 
forerunner generated by the onset of the source, and of a monochromatic 
front which carries the oscillation frequency of the source. For a medium 
with Schrodinger-like dispersion the monochromatic front propagates with a 
velocity which is in agreement with the traversal time, and in the relativistic 
case the velocity of the fronts is limited by the velocity of light. For sources 
with a sharp onset, the forerunners are not attenuated and in magnitude far 
exceed the monochromatic front. In contrast, for sources which are frequency- 
band limited, the forerunners are also attenuated and become comparable to 
the monochromatic front: like in the propagating case, there exists a time at 
which a broad frequency forerunner is augmented by a monochromatic wave. 

PACS numbers: 05.45.+b, 72.10.Bg, 72.30.+q 
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Fifteen years ago Rolf Landauer in collaboration with one of the authors [|IJ revived an 
old question: What is the time of interaction of a tunneling particle with the barrier? This 
work was motivated by the insight that the most often used approach which follows peaks 
of wave packets as they approach and emerge from a tunneling barrier is dubious: there 
is no causal relationship between the peak of an incident wave packet and the peak of an 
emerging wave packet. Instead a novel approach was needed: Ref. 1 and subsequent work 
investigated tunneling through a barrier with an oscillating potential and analyzed 
tunneling as a function of the oscillation frequency. Similarly, the Larmor clock, originally 
proposed by Baz || and Rybachenko ]5[], based on the precession of the spin of particles was 
re-analyzed and the importance of spin rotation was pointed out ||. Both the oscillating 
barrier and the Larmor clock lead to traversal times which differ from a stationary phase 
analysis. We can indicate here the wide interest and the broad discussion which these works 
have generated only with a reference to a few recent works and by a number of reviews 



While there have been only a few experiments on electrical systems |12| , recently a num- 



ber of interesting experiments which exploit the analogy ||T3| between electron tunneling 
and tunneling of photons |14| or the propagation of classical microwaves in evanescent me- 
dia |T5|-|nj] have provided a new impetus to the field. Of particular importance are the 
apparently superluminal tunneling times reported in these experiments. The superluminal 
appearance of wave packet peaks can be explained according to Deutch and Low |20| by 
demonstrating that the transmitted wave packet is made up mainly of contributions which 
stem from the head of the incident wave packet. Japha and Kurizki |2]J explain that the 
connection between the transmitted wave packet and the leading portion of the incident 
wave is a consequence of destructive interference of causal retarded tunneling paths. It is 
clear that the wave evolution is inherently a causal process and that superluminal velocities 
arise from the use of an acausal definition of the velocity like the comparison of the peaks 
of the incident and the transmitted wave. These experiments thus raise the question of 
whether it is possible to define a traversal time which represents a manifestly causal process. 
We can identify a number of different velocities which characterize wave propagation like 
the phase velocity, the group velocity, the signal velocity and the velocity of the head of the 
wave. It is well understood that both the phase velocity and the group velocity can exceed 
the velocity of light, but that the signal velocity is always smaller than or at most equal 
to the velocity of light and the velocity of the head of the wave always coincides with the 
velocity of light. In principle there should exist a traversal time corresponding to each of 
these possible definitions of velocity. It is the purpose of this work to discuss a traversal 
time which in an evanescent medium corresponds to a signal velocity, and which, therefore, 
characterizes a causal process. 

Interest in the signal velocity of propagating waves arose from the observation by R. Wien 
that, in apparent contradiction to special relativity theory, the phase velocity and the group 
velocity can easily exceed the velocity of light in dispersive media. In response A. Sommerfeld 
started a fundamental discussion by analyzing the propagation of wave fronts. While the 
phase velocity and the group velocity can exceed the velocity of light without contradicting 
the principles of causality and relativity, the velocity of the front and the signal velocity 



remain strictly smaller than or equal to the velocity of light [23]. A detailed discussion of 



this approach was subsequently given by Brillouin pi. The method of Sommerfeld and 



2 



Brillouin considers a source which is quiescent up to a given instant, when it is abruptly 
switched on. In a medium that allows propagation, the sudden onset of the source leads 
at the observation point located some distance away from the source to small forerunners 
of indefinite frequency; it is followed by a wave which oscillates with the frequency of the 
source and marks the arrival of the signal. For freely propagating waves (classically allowed 
regions of dispersion) the forerunners are small and the arrival of the main part is marked 
by a rapid increase in the intensity of the wave. The method of Sommerfeld and Brillouin 
was applied to evanescent media with a Schrodinger-like dispersion by Stevens [2^] . Stevens 
calculates an arrival time of a signal. He did, however, not analyze the magnitude of the 



different contributions to the total wave. Indeed, subsequent works by Teranishi et al. [26 



Jauho and Jonson []27j, Ranfagni et al. |28j have called into question the very existence of 



a front which would mark the arrival of the signal in an evanescent medium. Moretti [29 
emphasized exact analytical solutions, but similarly to Stevens did not discuss the magnitude 
of different contributions to the evolution of the wave. The lack of a main part found from 
numerical analysis p7| , |2~8|| was further emphasized in more recent analytical work by Brouard 
and Muga |5C|| . In contrast to Stevens who predicted an arrival time of the main part of 
the wave with a traversal time which is in agreement with that found from the oscillating 



barrier or the Larmor clock, Ranfagni et al [2"S] and Brouard and Muga BO] suggest that 



their analysis is actually in better agreement with the phase time. Since propagation occurs 
now in an evanescent medium, it is essential to investigate the amplitudes of the front and 
the forerunners in detail. Contrary to the case of freely propagating waves, the evanescent 
medium transmits the high-frequency components which make up the forerunners with little 
attenuation (the forerunners propagate in effect freely) while the "main part" of the wave 
is evanescent and thus exponentially suppressed. 

The situation changes dramatically if instead of a source with a sharp onset a frequency- 
band limited source is used. The source which is of interest switches on gradually but still 
fast compared to the traversal time. For a source which switches on too slowly, the traversal 
time cannot be resolved. Moreover, the highest frequency must be smaller than the threshold 
which permits free propagation. Thus the frequency-band is limited both at the high end 
and at the low end. For such a frequency-band limited source, all frequencies are in the 
evanescent range of the dispersion. This has the consequence that not only the fronts which 
carry the monochromatic frequency of the source, but also the forerunners are exponentially 
attenuated! The forerunners may, however, still exceed in magnitude the monochromatic 
fronts since the evanescent waves of the highest frequencies of the source are attenuated less 
strongly then the evanescent wave with the frequency of the source. Within this limitation, 
we encounter for the frequency-band limited source a wave evolution which is completely 
analogous to the freely propagating case: at the traversal time known from the oscillating 
barrier |l| and the Larmor clock ||, a forerunner of indefinite frequency is augmented by 
a nearly monochromatic wave of comparable amplitude. The discussion for the frequency- 
band limited source is presented for the case of a particle with Schrodinger-like dispersion. 
The discussion for a relativistic dispersion will be presented elsewhere 
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Source with a sharp onset 



We want to investigate propagation of a particle field ^(x,t) into a region x > with 
a constant potential V. The field equation — Schrodinger equation for a non-relativistic 
particle, Klein-Gordon equation for a relativistic boson, Dirac equation for a relativistic 
fermion — has plane-wave solutions 

^(x,t) = e- l[ujt ~ lk{uj)x] (ureal), (1) 

where the wave number k(u) is determined by the dispersion relation 

h 2 k 2 

Hu = V H (non-relativistic particle), (2) 



(tiu — V) = (mc ) + c h k (relativistic particle). (3) 

k(u) is real for propagating waves and imaginary for evanescent waves. Its sign is determined 
by the boundary condition that il> u (x,t) is an outgoing wave in the +x-direction: Positive 
group velocity for propagating waves and exponential decay with +x for evanescent waves 
yields 

k(u) sign(o;) > for propagating waves (4) 
Imk(u) > for evanescent waves. (5) 

For given u, the dispersion relations (QJ^) have only a single root satisfying the boundary 
condition Eq. (flJ^). 

Following Sommerfeld and Brillouin |23 , 24]] and Stevens [25|] we consider an arrangement 
which permits the investigation of wave front propagation: For t < 0, the field is everywhere 
zero, 

tp(x,t) = (t<0). (6) 
At t — 0, a source located at x = of frequency Uq and amplitude A is switched on, i.e., 

V>(0, t) = A{t) e~ iuJot with A(t) = A 0(t), (7) 
where Q(t) is the step function, or in o;-space, 

^(0, u) = A(u - u ) = / A(t) dt = — — . (8) 

Jo u — ujq + «U + 

The solution satisfying these initial conditions is 

$(x,u) = A{u - u ) e ik(uj)x (x>0), (9) 

or in t-space, 

1 p-\-oo 

^(x, t) = — A(u- u ) e -<M-*H»l du. (10) 
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For the following, it is convenient to introduce a frequency corresponding to the kinetic 
energy (plus rest energy in the relativistic case) of the particle, 



hfl = hu-V. (11) 
Then, the solution may be written in the form 

Mx, t) = — e~* vt / — — -± e-W*.*) dH, (12) 

^ v ' ; 2vr J-oo n-n + io + K ' 

where the phase function (f>(ft; x, t) is given by 

(p(n;x,t) =Qt - k(Q)x. (13) 

The integration is along the real fi-axis above the pole at Q = Q$ — + and above the branch 
cuts of the phase function 0(f2; x, t). The strategy consists in shifting the path of integration 
into regions of the complex f2-plane where Im0(f2; x, t) — > — oo as \Q\ — ► oo and collecting 
any contributions from the pole and the branch cuts. 



Non-relativistic particle 

The phase function 4>(Q) in the complex f2-plane is shown in Fig. [I]. It has a branch cut 
Q = (0 . . . oo). The values of the function </>(f2) on the lower sheet are related to those on 
the upper sheet by 



flower ^ = [^«PP«r(n*)]*. 



(14) 



The phase factor e has a saddle point on the real axis. The saddle-point condition 
d(j)/d£l\ s = yields 



HQ* 



m x 
~2¥ 



x 



hk s — m — , h(p s = —hQ s t = — 



m x 



t ' " 2 t 

and the second derivative of the phase function at the saddle point is 



d 2 



dQ? 



t 



2a 



h t 3 



m x A 



(15) 



(16) 



The line stph of stationary phase Re0(f2) = <p s , on which Im0 decreases to — oo as one 
moves away from the saddle point, is the parabola 



a. 



(17) 



where we have introduced the notation Ref2 = Q r , Imfl = On the stph-line, the phase 
function is given by 
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FIG. 1. Form of the phase function <j)(£l) in the complex $7-plane for a non-relativistic particle. 
Shown are the lines of constant <p r = ~Re<J)/(f> s for ip r = —4... + 6 (underlined numbers) and 
ifi = Iva.(j)/(j) s for ifi = —6 . . . + 1. The branch cut is indicated by an increased line thickness. The 
part of the line of stationary phase Re <j) = <p s on the upper and lower sheet is shown as a full and 
a dashed thick line, respectively. After Ref. [[H]] . 



1 + 



1 - 



a 



(18) 



For a fixed point x in space, the crossing points Q = ±fi s of the stph-line with the real axis 
move inwards with increasing t. One of them crosses the pole Q = Qq at time t = x/v m 
where velocity given by p5| 



— h\n \ 



m 



(19) 



We now deform the path of integration away from the real f2-axis such that it coincides 
with the stph-line. As long as t < x/v m there is no obstacle, and the integral remains 
unchanged; but for t > x/v m there occurs a loop around the pole which gives rise to a 
contribution ip p to the wave function of the form 



tP P (x, t)=A e-'^+^y e lk(n ° )x e(v m t - x). 
The integral along the stph-line yields a contribution tp s to the wave function, 



ip s {x,t) 



iA 
2^ 



J s 



1 



tph fl — f2( 



_ 1 

e 2 



dn. 



(20) 



(21) 



It may be evaluated in Gauss approximation, 
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if the width ((f)'^)" 1 ^ 2 of the exponential phase factor at the saddle point is small compared 
to the distance over which the prefactor changes significantly, 

0l' (n. - fi ) 2 = I (^^-) 2 n s t » 1. (23) 



2 v n s 

The pole contribution (|20|) describes a monochromatic wave with the frequency uo = u of 
the source, with a front which travels with velocity v m given by Eq. (|T§|). In the propagating 
case frfi > 0, v m is equal to the group velocity v g = du/dk; in the evanescent case hfl < 0, 
v m is related to the traversal time for tunnelling Q r for large barrier widths L by v m = L/r. 

The saddle contribution fl2~ID , on the other hand, describes the perturbation due to the 
switching-on of the source at t — which contains arbitrarily high frequencies; its front 
travels with the maximum speed permitted by the wave equation, which is infinite for a 
non-relativistic particle. 

Pole contribution and saddle contribution correspond to the "main part" and the "fore- 
runner" of Sommerfeld and Brillouin, respectively [p3p4f . We will keep the expression 
"forerunner" for the saddle contribution, but will call the pole contribution the "monochro- 
matic part" instead of "main part" , because the latter expression is extremely misleading in 
the evanescent case where the monochromatic part is exponentially small compared to the 
forerunner. 

It is important to point out that the existence of a wave front of the monochromatic 
part ip p does not mean that the total wave function is discontinuous at x — v m t. In fact, 
the solution of the Schrodinger equation is continuous for all x > 0,t > 0. Thus, the pole 
contribution ip p must combine at x — v m t with the saddle contribution ip s in such a way 
that the total wave function ip is continuous. 

This continuity requirement leads to an independent determination of the front velocity 
v m of the monochromatic part: The condition that the real parts of the phases of ip p and i/j s 
coincide for x = v m t reads in the propagating case (S7o > 0, k(Qo) real) 

2 

Tfl X 

Q t - k(Q Q )x = - — ^ for x = v m t, (24) 

and in the evanescent case (Qq < 0, k(flo) imaginary) 

2 

TTX X 

- \Q \t = -w±-r fo r x = v mt, (25) 



which yields the same result as Eq. fll9|). This is important because one can raise an objection 
against the derivation based on the crossing of the stph-line with the pole: Although the 
stph-line is the natural choice for the saddle-point integration, this choice is by no means 
unique: the integral in Eq. ( |2T|) does not change if the path of integration is shifted away 
from the stph-line. 

We now show explicitly that the crossing of the stph-line with the pole gives rise to a 
discontinuity of the saddle contribution i/; s (x,t) at x = v m t which exactly compensates the 
discontinuity of ip p (x,t) at the onset of the monochromatic wave. 
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In the propagating case, in the integral of Eq. (|2ip we set Q r = Q s (l + u) where u is 
an integration variable which measures the distance from the saddle point Q s on the real 
f2-axis. ^From Eq. (|T7D, we obtain 

Q = Q s (l + (1 -i)u- l -u 2 ^j ; dn = Q s (l - i - iu) du. (26) 
Then, the integral I along the stph-Yme becomes 

J-oo (1 — l)U — — Uq 

where we have introduced the abbreviations 

Q — Q s (v m t) 2 — x 2 „ x 2 

where Ax = h/(mv m ) equals (up to a factor of 27r) the de Broglie-wave length of a par- 
ticle moving with velocity v m . The parameter in the validity condition ( ^3|) for the Gauss 
approximation takes the form 

n s -n \ 2 ^ 2 ( vjt 2 -x 2 \ 2 Vm t 

-nT) nst = au ° = {—^-r ) 2Kx- (29) 

The Gauss approximation fl22|) remains valid as long as au 2 ^> 1, which is satisfied except 
very close to x — v m t. In the immediate neighborhood x = v m t where au^ < 1, on the other 
hand, the integral becomes 

r+ui \ 

1= lim / t-. r — du = iir signw , (30) 

«i-oo7_ Ul M _ 1(1 +l ) UQ 6 V I 

which yields the saddle contribution 



-Ae R V 2 « J s ign( 



^ s (x,t) = --Ae H 2 « ; sign(w m t - x) (v m t -> x). (31) 

Thus, in the propagating case, the saddle contribution has at x — v m t a jump of magnitude 
A. 

In the evanescent case, the Gauss approximation (|22| ) remains valid at x = t> m £ where 
fi — —Q s - However, since the wave number k(Q ) is imaginary, the "main part" ip p is 
exponentially small compared to the "forerunner" ip s , and must therefore be compared with 
the equally small contribution to the integral in Eq. ([H]) from a narrow region around the 
pole. We set fl r = Q s (—l + it) where the integration variable u measures the distance from 
the point — Q s on the real f2-axis, and Qq = — O s (l + wq) such that 

l^ol-a vlt 2 -x 2 tt -tts v^t 2 + x 2 
w = ^ = - 2 whereas u = — = . (32) 

^From Eq. (fPT|), we obtain 
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tt = tt s i^-l + (I + i)u - \ 2 J ; dQ = Q s (l + i - iu) du. (33) 

The contribution to the integral from the narrow region —U\ < u < Ux, where \wq\ <C Mi <C 
1, au\ <C 1, is given by 

AI = e~ 2a [ +U1 y— r — du = iir e~ 2a sign w (34) 

J-ui u + ^(1 — 

which yields the contribution to ip s 



I — — I vt— — — 1 ^ 

Ai/j s (x,t) = — Ae H 2 W e "^ a Tsign( 
2 



) m t-x) (v m t^>x). (35) 



Thus, in the evanescent case, the jump of the saddle contribution at x = v m t is exponentially 
small. 

Comparison with Eq. ( |2DD shows that ip s has indeed in both discontinuity opposite 

to that of ip p . 

To summarize: In the propagating case the wave as observed at a distance x from the 
source grows to one half of the asymptotic amplitude A at a time t = x/v m when the 
monochromatic (main) contribution of the source sets in. In contrast, in the evanescent case 
the sharp onset of the source generates a contribution to the wave which largely determines 
what is seen at the observation point. The saddle point solution is of order A at t = x/v m 
when the exponentially small monochromatic contribution arrives. For a Schrodinger-like 
dispersion there is no limiting velocity, the sharp onset of the wave generates immediately 
a wave at the observation point x. Next we investigate a medium with dispersion which 
contains a limiting velocity. 



Relativistic particle 

The phase function 0(f2) in the complex f2-plane is shown in Figs. |^ and |3| for x > ct 
and x < ct, respectively. It has branch cuts hQ = (—00 . . . — mc 2 ) and hQ = (+mc 2 . . . +00), 
and the values of 0(fi) on the upper and lower sheet are again related by Eq. ([141). 

For x > ct, Im0(f2) — > —00 for ImQ — > +00 (see Fig. H). By shifting the path of 
integration upwards, one recognizes that 

ip(x,t)=0 for all (x, t) with x > ct, (36) 

as required by relativistic causality. 

For x < ct, on the other hand, lm<f)(Q) — > —00 only for Imfl — > —00 (see Fig. |^). In this 
case, the phase factor e~^ n ^ has saddle points = ±fi s on the real axis. For the saddle at 
+£l s , the saddle-point condition d<p/dQ\ s = yields 

t x $ 2 

hQ s = mc 2 —, hk s — m — , h(p s = hQ s — = mc 2, d, (37) 

where we have introduced the quantity 
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X 2 



d(x,t) = \t 2 --. (38) 



c 2 



The second derivative of the phase function at the saddle +Q S is 



h $ 3 

(39) 



m x 2 



The line stph + of stationary phase Re0(fi) = <p s through the saddle +Q S , on which Im0 
decreases to — oo as one moves away from the saddle point, is given by 

Q = (Q r -Q s )[h 2 Q r Q s -(mc 2 ) 2 ] 

^[h 2 n 2 - (mc 2 ) 2 ]{-h 2 n 2 + 2h 2 n r n s - (mc 2 ) 2 ] ' 

It crosses the real axis at the saddle point Q = Q s and at Q = (mc 2 /h) 2 /Q s , and goes to 



f2j = — oo for Q r = il s ± yf2 2 — [mc 2 /K) 2 . In contrast to the non-relativistic case, it is no 
longer possible to obtain the imaginary part of the phase function <p(Vl) on the stph-Ym.es 
in closed analytic form. At the crossing point Q = (mc 2 jK) 2 /Q s , its imaginary part has the 
value lm<f) = —imx 2 /ht, which is the same as the value of Im0 in the nonrelativistic case 
at the crossing point Q = —Q s . 

The results for the saddle at — Q s are obtained from Eqs. (|3?| - |40"1) by substituting Q r \— > 

For a fixed point x in space, the crossing points of the stph + -\ine with the real axis move 
towards Q = mc 2 as t increases form x/c to oo. In the propagating case, the saddle point 
Q s crosses the pole Q = Qq > mc 2 /h at time t = x/v m where v m is given by 



mc* 



\ \nn 



(hn > mc 2 ), (41) 



which agrees with the group velocity v g = du/dk. In the evanescent case, the point Q 
(mc 2 /h) 2 /(Q s ) crosses the pole Q = Qq < mc 2 /h at time t = x/v m where v m is given by 



\ 




(hn < mc 2 ), (42) 



which is equal to h\k\/m. We relate v m to the relativistic traversal time for tunneling Tt r for 
large barrier widths L by v m = L/r. 

The velocity v m exceeds the velocity of light c neither in the propagating nor in the 
evanescent case; it becomes equal to c only for zero energy and in the limit of infinite 
energy The energy dependence of v m is shown in Fig. f|. 

In order to identify the contributions to the wave functions tp(x,t), we proceed in the 
same way as in the non-relativistic case. The path of integration is deformed away from the 
real il-axis such that it coincides with the stp/z ± -lines at positive and negative energies. As 
long as t < x/v m , no obstacle is encountered, and the integral remains unchanged; but for 
t > x/v m there occurs a loop around the pole which gives rise to a pole contribution 
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Re Q/Q s 

FIG. 2. Form of the phase function (p(Q) in the complex f2-plane for a relativistic 
particle for ct = 0.75x. Shown are the lines of constant (p r = K Re cj) / (mcV x 2 — c 2 t 2 ) for 
cp r = —6 . . .+6 (Oj < 0), =0, ±0.5 (Oj > 0) (underlined numbers) and ipi = % Im 4>/ (mcy/x 2 — c 2 t 2 ) 
for W = -4, -2, -1.5,-1, -0.5. After Ref. EB. 






Re Q/Q s 

FIG. 3. Form of the phase function 0(0) in the complex $7-plane for a relativistic particle for 
ct = 1.25x. Shown are the lines of constant ip r = hUe (f> j \mc\J 'c 2 t 2 — x 2 ) for (p r = —6 . . . + 6, and 
ifi = h Im (f>/(mcV c 2 t 2 — x 2 ) for ipi = —6. . . + 1. The branch cuts are indicated by an increased 
line thickness. The parts of the lines of stationary phase Re^ = 4>f on the upper and lower sheet 



are shown as full and dashed thick lines, respectively. After Ref. [31| 
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0.0 — 




(/zco -V)/mc 2 - 

FIG. 4. Front velocity of the monochromatic wave for a relativistic particle as function of 
energy. After Ref. Q. 

ij p (x,t) = Ae-<* v+Qo > e lk{Uo)x &(v m t - x) } (43) 



describing a monochromatic wave with front velocity v m given by Eqs. (|41|,42). The integral 
along the stph^-lines yield saddle contributions 

tfUA = ile-'U™^) / -l—e-^ n UQQ(ct-x). (44) 

Z7T Jstph ± 11 — \Iq 

with a front which moves with velocity c. The contribution due to the saddle at +Q S describes 
the excitation of particles, and that due to the saddle at — fl s the excitation of antiparticles. 
The latter part will be significant only for energies deep in the classically forbidden region, 
HQq ^ ttlc 2 . The saddle contributions may be evaluated in Gauss approximation, 



tf{x,t) = ^ tF^-fr TT^TT e"^ i±mc ^ Q{<* - x). (45) 

Ys v ' 1 2vr V mc 2 i? c 2 t 2 fi s ± fi v ; v ; 

if 

4' a (fi, =F fi ) 2 = (Q s =F fi ) 2 » 1. (46) 

Since the solution of the wave equation is continuous for all x < ct, the pole contribution 
i/jp must combine at x — v m t with the saddle contribution ip s in such a way that the total wave 
function ip is continuous. This requirement leads again to an independent determination of 
the front velocity v m of the monochromatic part: The condition that the real parts of the 
phases of ip p and ip s coincide for x = v m t reads in the propagating case (hflo > mc 2 , k(Q ) 
real) 

h (Q t — k(Q )x) = mcV c 2 t 2 — x 2 for x = v m t, (47) 
and in the evanescent case (M7 < mc 2 , k(Q ) imaginary) 
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hil t = mc\/c 2 t 2 — x 2 for x = v m t, 



(48) 



which yields the same result as Eqs. PT| . [T^| ). 

Like for a non-relativistic particle, it can be shown explicitly that the crossing of the 
stph + -line with the pole gives rise to a discontinuity of the saddle contribution ip+(x,t) 
at x = v m t which exactly compensates the discontinuity of i[) p (x,t) at the onset of the 
monochromatic wave. One obtains in the propagating case 

^f(x,t) = ~Ae~^ vt+mc2 ^ signet - x) (v m t -> x). (49) 

and in the evanescent case 

A^+(x,t) = - l -Ae-^ vt+mc2 ^e-^ signet - x) (v m t -> x). (50) 

Like for the non-relativistic particle, in the evanescent case the saddle point contribution 
far exceeds the monochromatic contributions to the total wave. 



Fronts of a frequency-band limited source 

Clearly, the high frequencies generated by a source which is switched on instantly are 
highly undesirable. In this section we investigate a source which is limited in the frequency 
band width. The fact that experimental signals are often frequency-band limited has already 
been emphasized p2|j33| (that is of course a technical and not a fundamental limitation 



We restrict the discussion to the case of a non-relativistic particle. The frequency- 
band limited source is described by an amplitude in Fourier space which is the product 
A(uj — ujq)x(uj) of the amplitude A = iA/(u — ujq + i0 + ) of the source with a mathematically 
sharp onset and a characteristic function x(oS) which limits the range of frequencies. 
For simplicity, we take 

X(oo) = &(uj - (uj - Aoj)) - Q(oj - (uj + Aw)), (51) 

where is the step function, although strictly speaking this form violates causality: As can 
be seen from the Fourier transform, 

, (f) = i^.^M, (52) 

7T t 

the source now radiates already for t < 0. The frequency band width Au is chosen such 
that 

1 V 

- < Aa; < — -ujq = \tt Q \. (53) 



T 



h 



Here the lower limit ensures that the onset of the source is still fast compared to the traversal 
time t = x/v m and that the effect of causality violation is kept mall. The upper limit assures 
that all frequencies contained in the source are in the evanescent region. 
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We start with the integral for the wave function given by Eq. (|12|). The integration path 
is originally along the real axis from = Q — Au to Q + = Qq + Au; above the pole at CIq. 
It is now convenient to deform the integration path in such a way that the path first follows 
the line of constant imaginary phase Im</>(0_) = $\ which goes through the point then 
from the point where this line (see Fig. [I]) intersects the stp/t-line follows the line stph to 
the point where it is intersected by the line of constant imaginary phase Im0(f2 + ) = <pf, 
and then follows this line to the point + . For Qq < —Q s , we have to pull the integration 
path across the pole, which yields a monochromatic contribution to the wave function 

%{x, t)=A e -<s v + n o)t e ik(n )x Q{Vmt _ x ^ (54) 

of the same form as Eq. ( p0|) . 

The other three contributions to the integral are evaluated for a time interval for which 
—Q s is close to the pole Qq. Furthermore, we are interested in sources with a band width 
which is small compared to |Oo | . To investigate this time window, we can therefore linearize 
the phase function around the frequency — Q s . In terms of the variables u and v which 
measure the deviation from Q s , i.e. Q = Q s (—1 + u + iv), we find for the phase function 
<f) r = Re0 = Q s t(—1 + u — v) and 0j = Im0 = Q s t(— 2 + u + v). Locally, around Q = —Q s 
the lines of stationary phase form thus a rectangular grid. In terms of u the upper and lower 
frequency limits are 

Q s - \Q \ Au Auj 

"± = ^1^ 717 = -* ±7)7- (55) 

We recall from Eq (^) that Wq = {v^t 2 — x 2 )/x 2 . The first part 1\ of the integral along the 
line 0j~ = Q s t(—2 + uJ) extends in u from U- to u_/2, the second part I2 of the integral 
along the stph-line with <p r = —Q s t extends in u from w_/2 to u + /2, and the third integral 
^3 along <pf = Q s t(—2 + u + ) extends in u from u + /2 to u + . 

The first (index — ) and the third integral (index +) is approximated by pulling the 
denominator for u — u± in front, 

MXit ) « (±±M J_ e -*K^) e -f ^ e H+Wstu ± sy n (n st ^±) . (56 ) 
2ti Au t \ 2 J 

For | | 3> Au;/|r2o|, which holds outside a narrow time interval around t — r, the second 
integral is approximated by pulling the denominator at u — —wq in front, 



'1 + i)A 1 -j-( vt "' 

: : — e 1 V 



ip 2 ( x ,t) * - — e H'" 2 ?)e-tte^-W*swh{Au,t). (57) 

Z7T (lis — I i Zq I J ^ 

The key point is that all three integrals are exponentially suppressed. For t = r we have 
exp(— mx 2 /{2Tir)) = exp(—mv m x/h). The second point to notice is that because we still 
have a range of frequencies, the uppermost frequencies are least suppressed. Thus at t — r 
the upper frequencies are enhanced by a factor exp(Q s ru + ) = exp((mv m x/h)Au/\Q \) 
whereas the frequencies at the lower end of the spectrum are additionally suppressed by 
a factor exp(fi s r-u_) = exp(— (mv m x/h)Au/\Qo\). Clearly the different exponential sup- 
pression of these frequencies is unavoidable. 
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Let us next discuss what happens in the narrow time interval around t = r where 
| Wo | <C Ao>/|fio|, and where the expression for the wave function ip 2 (x,t) obtained above is 
not valid. For this integral, we have to reconsider a discussion analogous to that which leads 
to Eq. (34). For w ^ 0, by pulling the exponential factor at the upper limit in front, we 
find 

ip S 2(x,t) m -e S V 2 ' J e"--e Awt arctan — . (58) 

71 \U S - \il \ J 

At Wo = 0, on the other hand, the integral has a discontinuity which exactly compensates 
the jump of the pole contribution ip p . 

The discussion presented here does not show explicitly the short-time and the long- 
time behavior of the stationary phase solution. This restriction is due to the linearization 
of the phase function. If we are sufficiently far away from the point Qo = —Q s , we can 
estimate the contribution from the saddle-point solution near the intersection point of the 
stph-line with the vertical line that runs through Q + . This leads to a wave function whose 
magnitude is governed by exp(— ~Q s t(l — (Q + /Q s ) 2 )). For times which are short compared 
to r, the saddle-point frequency Q s is large compared to and the exponential can be 

approximated by exp(— |fl 8 t). For times which are large compared to r, we have Q s <C \ 
and the exponential can be approximated by exp(— ^Q s t(Q + /Q s ) 2 ) which decays with time 
as exp(— const, t 3 ). Thus the saddle point contribution is small both at times which are 
short compared to r and at times which are long compared to r. 

To summarize: By limiting the frequency band width of the source we obtain a wave evo- 
lution in the evanescent medium which shares the essential properties of the wave evolution 
in the propagating case: At a time t = r a forerunner with a broad frequency distribution 
is augmented by a monochromatic front with the frequency of the source. 

Discussion 

In this work we have focused on the question of whether it is possible to find a traversal 
time which is associated with a causal process. The process investigated is the propagation 
of monochromatic fronts. For the case of a source with a sharp onset these monochromatic 
fronts are exponentially small compared to the forerunner which is not attenuated. How- 
ever, as we have shown, if the source is frequency-band limited the forerunners are also 
exponentially attenuated. For the evanescent case, a wave evolution is found, which as in 
the propagating case, exhibits at t = r a crossover from a wave with a broad frequency 
spectrum (forerunner) to a nearly monochromatic wave. We also emphasize that the total 
wave function, respectively, the resulting probability distribution is not the only observable. 
Possibly, an other way to make the monochromatic fronts visible is to incorporate a detec- 
tor which is sensitive only to a narrow range of frequencies in an interval around the main 
frequency ujq of the source. Or, we could consider particles with a spin in an evanescent 
medium with a weak magnetic field and could investigate the rotation of the spin as a func- 
tion of time similarly to the stationary analysis of the Larmor clock M. Our work shows 
that the investigation of the time evolution of particle fields is a very intersting avenue of 
research. The comparison of such investigations with the results from the oscillating barrier 
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approach and quantum clocks [|l],|||8| should be particularly illustrative. Furthermore, we 
hope that this work stimulates also experiments which aim to identify the monochromatic 
fronts discussed in this work. 
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